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Abstract 

Recent observations with Swift have begun to uncover y-ray transients whose total energies are comparable to those of gamma-ray 
bursts (GRB), but have a duration an order of magnitude or more longer than the bulk of the GRB population. Some are suggested 
to form a new population of ultra-long GRBs, with a mean duration around 10 4 s, while a further population with y-ray durations 
> 10 5 s may represent manifestations of relativistic outflows from stars shredded around massive black holes in tidal disruption 
flares (TDFs). Here I review the observations of these new classes of events, discuss progress towards identifying their progenitors 
and suggest how new observations may both hone our understanding of the outbursts, and allow them to be used as probes, that 
offer both complementary and additional tools to GRBs. 

Keywords: Galaxies: Active; Gamma-ray bursts: general; Gamma-ray burst: individual (GRB 101225A, 111209A, 121027A, 
130925A); supemovae:general 


1. Introduction 


Many of the most energetic high-energy astrophysical sources 
show marked variability on a range of timescales. From short 
duration gamma-ray bursts (GRBs) to outbursts from active galax¬ 
ies (AGN) it is clear than the high-energy sky varies on timescales 
from milliseconds to many years. The extremes of this regime 
have been probed since the beginnings of X-ray and y-ray as¬ 


tronomy; GRBs were first recognised in the late 1960’s (Klebe- 


sadel et ah] |1973| ), and AGN were amongst the first identified 
extragalactic X-ray sources ( Elvis et ah) 1978| ). However, prob¬ 
ing durations of hours to days has proved to be challenging. 
Such systems are too rare and fleeting to be picked up by nar¬ 
row field X-ray telescopes (which are often used for studies 
of longer lived events) while they are also often too faint to 
trigger wide field rate-based y-ray detectors, which have been 
the dominant route for mapping the rarest, but brightest events 
(such as GRBs). 


The Swift satellite ( [Gehrels et al.[(200 4| offers an ideal route 
to the study of transient events in these time frames. Its wide 
field covers roughly 1/6th of the sky in a single view of its Burst 
Alert Telescope (BAT). It can be triggered not only based on the 
rate of incident photons, but by reconstructing an image based 
on a long integration, and searching this for sources that would 
not have triggered a rate-based detector. Indeed, it can pro¬ 
duce integrated triggers over timescales of several days, and so 
long duration transients have been identified both as GRB trig¬ 
gers, and via an ongoing transient monitor ( jKrimm et al. 2013 ). 
This is coupled with the rapid X-ray and UV/optical follow-up, 
meaning that these events, once found, can be localised and de¬ 
tailed multi wavelength follow-up initiated, both from ground 
and space-based facilities. 

Remarkably, it has only been in the last few years (the sec¬ 
ond half of Swiff s lifetime, at the time of writing) that such 


populations have begun to be uncovered, perhaps in part be¬ 
cause events in the early years may have gone otherwise un¬ 
recognised, but also because they are likely intrinsically rare. 
In this review I outline Swift's progress in unveiling new pop¬ 
ulations of transients that emit y-rays over time periods from 
hours to days. These duration are orders of magnitude longer 
than typical GRBs, but still much shorter than AGN-outbursts. 
I concentrate on the most energetic events, whose total energy 
releases exceed ~ 10 53 ergs. These events split broadly into 
two categories, with durations of ~ 10 4 and > 10 5 seconds (see 
Figure [I]), they may be related, but perhaps more likely repre¬ 
sent a diverse set of events with significantly different progeni¬ 
tors. For each class of transient I will discuss the observational 
characteristics as a class, and proceed to consider the plausible 
progenitor models. There is significant overlap between their 
proposed progenitor systems, which break-down broadly into 
massive star collapses (of stars with larger radius than those 
creating long-GRBs), and tidal disruption flares. This in turn 
suggests that while the progenitors of each class may be dis¬ 
tinct, it is possible that both classes of progenitor are present in 
the observed populations. 


2. Determining a duration 

The duration of an outburst is a fundamental parameter in 
its description. In addition to providing a useful tool in the 
classification of a given outburst (e.g the observed distinction 
between short- and long duration GRBs (SGRBs and LGRBs) 
( [Mazets et al.| |1982[ |Kouveliotou et al.] |1993| )), the duration 
of a transient also provides hints on the activity in its central 
engine ( Zhang et al., 2006; Bromberg et al., 2013), and has been 
suggested to provide direct constraints on the size of long-GRB 
progenitors (Bromberg et al. 2012| ). 
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Figure 1: Parameter space for high energy transient phenomena, in the Tgo 
versus average luminosity place, from |Levan et al.| ( p014a) . Large numbers of 
long and short duration gamma-ray bursts can be seen, but only a handful of 
even longer events, lasting 10 4 - 10 6 s. These two classes, named ULGRBs 
and TDFs are the subject of this review. 


Traditionally the duration of a GRB is defined as 790 , the 
time period over which 90% of the total fluence is recorded. 
This is a simple and powerful definition, however it suffers from 
a lack of comparability between instruments which may have 
different effective passbands (and bursts may last for different 
durations at different energies) and from the sensitivity of the 
instrument. A more sensitive instrument will follow the decay 
of a burst for longer. Depending on its light curve this could in¬ 
clude significant additional energy, and hence increase the du¬ 
ration. Indeed, a comparison of the distribution of durations of 
bursts detected by Swift and CGRO/BATSE clearly shows that 
the softer response of Swift tends to omit the typically harder 
short bursts. There is also a suggestion that the long bursts it 
detects have larger Tg 0 values, again likely because of the softer 
response (Figure [2]). 

In addition to the distinction between long- and short-GRBs, 
another feature apparent in the distribution of GRB durations 
(Figure [2} is the appearance of a handful of bursts with dura¬ 
tions much longer than the majority of bursts observed. While 
very long bursts have been identified with durations of > 1000s 


by BATSE, BeppoSAX and Konus-WIND (e.g. Levan et al. 


2005[|Parshin et al.]|2008| ) they have not been studied in detail, 


and these events were generally still much shorter than some 
which have recently been detected by Swift that can be seen in 
gamma-rays for in excess of 10,000 s - this is the suggested 
population of ultra-long GRBs. The characteristics of this pop¬ 
ulation are yet to be strongly defined, but it is clear that there is 
a tail of “normal” bursts extending to ~ 1000s (see Figure [2]), 
and so durations of 5000s or more may be needed to lie clearly 
in the ultra-long population, we discuss this in more detail in 
section [3] A further population of outbursts has been identi¬ 
fied with durations of > 10 5 s, these have been seen either as 
repeated GRB triggers, or via the BAT transient monitor. The 
durations are so extreme (y-ray emission detected for > 10 6 s in 
some cases) that they clearly represent a separate population of 
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Figure 2: Left: The duration distribution of GRBs recorded by BATSE (black) 
and Swift-BAT (red). Pronounced differences can be observed at the short end, 
where BATSE records many more short bursts due to its harder response. There 
is also an apparent preference for Swift to record somewhat longer durations 
that BATSE. The population of ultra-long bursts are visible at the right hand 
side around 10 4 s. A log-normal fit to the Swift distribution from |Virgili et al.| 
{2013) is shown as the red dashed line. The ultra-long bursts lie well above the 
prediction of this model. Right: The distribution of T\, urst defined by [Zhang] 
|et al. (2014), along with a two log-normal fit to the data. Typically Tburst » 
Tgo, but it is striking that in addition to a central core that is well fit with a 
log-normal distribution there is a significant tail with prolonged central engine 
activity. 


bursts (and probably a separate progenitor) to the LGRBs, these 
very long events, candidate tidal disruption flares, are discussed 
in more detail in section [4] 


3. Ultralong gamma-ray bursts at 10 4 s 


3.1. A new class? Hints from duration distribution 

Ultra-long bursts may be difficult to discover since long- 
lived, but low flux bursts are readily missed. Indeed, the Swift- 
BAT normally only records the duration as measured immedi¬ 
ately following the trigger in so-called “event” mode. In the 
longest bursts the high-energy y-rays may continue to be de¬ 
tected in the following orbit, although it seems that with the 
Swift-BAI sample this is rare, with only ~ 15 events found in 
survey data surrounding the burst ( f Sakamoto 2014| . Even of 
these only a handful exceed durations of 5000 s, and most had 
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Figure 3: Two very different bursts highlighting the difficulty in accurately 
defining a duration. The upper panel shows GRB 050724, a burst that probably 
arises from the short GRB population and lies in an elliptical galaxy (SGRBs 
are notionally defined by 790 < 2 s, although in the case of GRB 050724 in ad¬ 
dition to a short spike there is lower-level extended emission lasting for ~ 100s). 
While the lower panel, shows GRB 111209A, a suggested member of the new 
ultra-long population of bursts. In each case different durations have been indi¬ 
cated. The classical is determined by the duration of the y -ray emission, in 
this case recorded by BAT immediately after the trigger, or in the case of GRB 
111209A recovered in subsequent orbits. This clearly misses some engine ac¬ 
tivity, and so may not be a good diagnostic. Alternatively, 7b ur st defines the 
total activity, but doesn’t distinguish between low and high levels, and can have 
very different values of 7%,burst (T 90 determined within 7b ur st)- 


been previously identified as candidate ultra-long events based 
on their prolonged, luminous X-ray emission. 

A crucial question is whether, on duration grounds alone, 
these bursts represent a new population of events to the tradi¬ 
tional GRBs. Different conclusions have been drawn on this. 
Virgili et aT7| ( |2Q13| ) argue that the overall distribution is consis¬ 
tent with underlying duration distribution of long-GRBs, which 
they find to be reasonably fit as a log-normal distribution with a 
mean of 1.47 (Tqo = 47s) and <x = 0.51. There is clearly limited 
statistical power at very long durations, where the small sample 
sizes result in significant counting errors. However, they note 
that the number of bursts observed with durations > 600 s (11), 
is consistent with expectations of the model (9.75). However, in 
this case we would expect the majority of the longest bursts to 
be observed closer to this cut (i.e. with durations of 600-1000 
s), in contrast, the ultra-long population suggested above has 
durations > 5000s. This lies > 4cr from the mean inferred by 


[Virgili etliL] ( |2013| , and given the number of bursts observed 
by Swift we would expect to observe none with Tg 0 > 5000s, 
compared with > 4 systems observed. This is shown graphi¬ 
cally in Figure [3] and may indeed suggest the bursts are a new 
population. 

Unsurprisingly, these difficulties in measuring durations have 
led to several authors investigating alternative tools for duration 
determination. Of particular importance is the discovery early 
in the Swift mission that the engines of GRBs are active for 
much longer than the prompt emission time, with a combina¬ 
tion of flares and plateaus indicating that in both long and short 
GRBs the engine can be active for many thousands of seconds 
beyond the end of the y-ray emissio n {jNouse k et al. 
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Zhan g et aL||2006| ). Motivated by this Zhang et al. (2014) have 


attempted to define an engine lifetime ( Tburst ) for many of the 
bursts detected by Swift. This is a difficult task since it requires 
some decomposition of the afterglow light and prompt emis¬ 
sion, both of which are likely to be contributing within a few 
hundred seconds of the burst. However, the cessation of engine 
activity is often followed by a rapid drop-off in brightness (for¬ 
mally more rapid than readily allowed by fireball blast wave 
models). By identifying this time it is possible to place some 
constraint on the engine activity time for 3/4 of the Swift bursts, 
although it is only relatively robust in ~ 1 / 2 of the bursts. The 
resulting distribution is shown in the lower panel of Figure [2j 
while the outline of the duration measurement method is shown 
graphically in Figure [2] It is notable in this distribution that 
the engine activity time is frequently much larger than Tgg. It 
is also clear that while with the exception of a handful of very 
long bursts the Tg 0 distribution is well fit by a log-normal distri¬ 
bution, the Tburst distribution consists of a much larger tail, that 


requires an additional component to provide a good fit (Zhang 


et al. 2014). Indeed, a similar approach, with a similar dura¬ 


tion distribution is reached by Boer et al. ( 2015), who define a 
duration, T x to be the last rapid decay in the X-ray lightcurve 
(a similar definition to Tburst)- However, [Zhang et al.| ( |2014| ) 
and |Boer et ah] ( |2015| ) show less agreement over the signifi¬ 
cance and interpretation of this tail, with Zhang et al. (2014) 
suggesting that an ultra-long population cannot be confirmed, 
while Boe r et aL\ (2015| ) infer that the ultra-long population is 
statistically different. Such diagnostics are complex, the long 
tail is not necessarily representative of a new population since 
there is no strong physical motivation for the log-normal dis¬ 
tribution of durations. Additionally, it should be noted that the 
rapid fall-off in the Tburst distribution occurs at around the Swift 
orbital period. Given that observations of bursts cannot nor¬ 
mally continue much beyond the first ~ 1000s it is expected 
that this could result in a pile up of Tburst values on either side 
of this gap ( [Zhang et al.||2014| . Nonetheless, that the distribu¬ 
tions of T 90 and Tburst are better explained with the presence of 
an additional population at larger durations, while being well 
described around their peak by log-normal distributions is sug¬ 
gestive of a new, additional population. To ascertain if these 
necessarily imply different progenitor systems it is necessary to 
turn to the multi wavelength afterglow properties. 
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3.2. The ultra-long GRB sample 


While there have been several examples of very long bursts 
over the years, the population as a whole remains poorly de¬ 
fined, and indeed is extremely difficult to define. However, 
the arguments above suggest that placing the duration cut to¬ 
wards the long end of the bulk of the GRB distribution (i.e. 
in the region 500-1000 s) will result in difficulties in isolating 
clean samples. Furthermore, it is the events which last much 
longer that are clearly distinct, and cannot be described with the 
same log-normal distribution that adequately describes long- 
GRB durations. Therefore, conservatively one could set the 
ultra-long bursts as those which have y-ray emission detectable 
at the sensitivity of current instruments (INTEGRAL, Fermi- 
GBM, Konus-WIND, BAT, MAXI etc) of ~ 5000s or mortQ 
This results in a small sample of bursts that can be considered, 
although it may be enlarged somewhat by the ongoing search 
for further events in BAT survey data ( [Sakamoto) |2014| ), or by 
the consideration of similar morphological properties in X-ray 
light curves that might indicate bursts belonging to the same 
population (Levan et al., |2014a] ). The events for which detailed 
mulitwavelength follow-up are available are shown in Table [I] 


They are GRBs 101225A, 111209A, 121027A and 130925A. 

However, while this route of identification is robust in the 
sense that it identifies bursts in which the prompt emission has 
lasted for a very long time, it may also miss other members of 
this class. The total fluence measured for the ULGRBs above 
is similar to that seen in most LGRBs, but spread over a longer 
time (~ 10 4 s versus ~ 10 2 s), this means that these bursts are 
harder to detect in y-ray s, and had they been fainter then they 
could easily have failed to have the longer lived y-ray emission 
detected, even if the brightest emission was sufficient to trig¬ 
ger the detector. Additionally, the orbit of Swift and Fermi is 
such that they cannot observed the position of the burst continu¬ 
ously for several thousand seconds. Indeed, Swift rapid slewing 
across the sky means that it spends a few hundred seconds on 
average at a given position. These satellites may then miss some 
of the later time y-ray emission. In the case of GRB 130925A 
the brightest episodes as recorded by MAXI or Konus-WIND 
happened at points where Swift was not observing the source, 
meaning that even if Swift data is considered over several orbits, 
it provides a poor view of the total duration of the burst. 

Therefore it is relevant to consider if other indicators may 
be used as diagnostics of the ULGRB population. These might 
be ( [Levan et ak|l2014a| ) 


• Prolonged X-ray emission (lasting ~ 10 4 s or more) within 
two orders of magnitude of the peak (i.e. at the start time 
of the X-ray observations). 


• A rapid decay at the end of this period, similar to those 
interpreted as high-lattitude emission in early GRB af- 
terlgows. 


'Note that by this definition, some bursts that have previously be ascribed as 
ultra-long, such as GRB 020410 (Levan et ah] 1 2005) and GRB 090124 (Virgili| 
|et al.|[20l3) are not in the sample. 


• Variability within the X-ray light curve through this pe¬ 
riod that is consistent with ongoing prompt emission, but 
unlikely to have an afterglow origin (e.g. rapid flaring, 
and a strong hardness-intensity correlation). 


Many more bursts meet one or more of the criteria above, 
although if these are truly representative of the ULGRB popu¬ 
lation is less clear. Indeed, the description above is very sim¬ 
ilar to that used by [Zhang et al.[ ( |2014| ) to define Tburst with 
the additional criteria that the luminosity of the burst must re¬ 
main within two orders of magnitude of the peak throughout 
this duration in order to classify as a ULGRB. A more quanta- 
tive approach would be to utilise Tburst as a measure of the total 
engine activity period (rather than use the total time over which 
y-rays are detected), and then determine T^ urst as the time 
over which 90% of the total fluence released during Tburst was 
recorded. In this case many GRBs would still have a short r 90 , 
even though Tburst was very long (e.g. GRB 050724, Figure 3), 
while because ULGRBs remain bright for longer they have a 
much larger value of T^^urst- Further investigations into this 
diagnostic may yield stronger constraints on the possible size 
of ULGRB population. 


3.3. Cosmological origin 

The properties of the first ULGRBs to be identified were 
sufficiently different from those of most GRBs that it was ini¬ 
tially unclear if they were of cosmological or Galactic origin, 
a debate very similar to that conducted for LGRBs and SGRBs 
in the early 1990s ( |Lamb| |1995| |Paczynski] |1995| ). For GRB 
101225A two very different models were put forward, the first 
that of a comet or asteroid being tidally disrupted around a neu¬ 
tron star ( [Campana et al.| |2011| ), the second the explosion of 
a supernovae inside a dense envelope (probably its own com¬ 
mon envelope, (Thone et al. [ |201 1| ). The former model indeed 
bares significant resemblance to a popular model for the cre¬ 
ation of Galactic GRBs ( Harwit & Salpeter j 1973), and has even 
recently been mooted as a possible origin for fast radio bursts 
(FRBs, |Geng & Huang ||20 15) , the latter model is a rare variant 
on binary channels to create LGRBs ( [Izzard etak| 2004[ Levan 


|et al.| 2006| ), in which the explosion occurs during the common 
envelope via the merger of a He core with a neutron star. In this 
model the afterglow evolution was reasonably fit with a type Ic 
supernova at z ~ 0.3 (Thone et al., |2011[ ) 

The redshift of GRB 101225A proved extremely challeng¬ 
ing to derive, since despite a bright afterglow (see below), it 
was largely featureless, showing no obvious emission or ab¬ 
sorption lines. However, deep, later observations taken with 
Gemini showed a series of faint emission lines at z = 0.85 
(Levan et al., [2014a) , resolving the Galactic versus cosmologi¬ 
cal debate, but placing it at a redshift well beyond the model of 
Thone et al.|( |201l| ). Subsequently redshift measurements have 


been rather more straightforward for GRB 111209A (z = 0.667 
from both absorption lines in the afterglow and emission lines 


from the host) ([Levan et al. 

2014a), GRB 121027A (z=1.71 

from afterglow absorption, (1 

fanvir et al. 2012), and host emis- 


the host galaxy emission lines ( [Sudilovsky et al. 2013| Q. Hence 
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the ULGRB population is clearly cosmological, and to date has 
a rather lower mean redshift (z ~ 0.9) than the LGRB popula¬ 
tion as a whole (z ~ 2.24 ( [Jakobsson et al.| [2006112012| )). This 
differing redshift is not surprising since the lower (on average) 
peak fluxes of ULGRB s would naturally result in a more limited 
horizon. 



Time since outburst (restframe days) 


Figure 4: X-ray luminosity light curves of various categories of X-ray emitting 
transient, including supernovae, long-GRBs and short-GRBs. The new classes 
of very long lived and energetic transient are shown at the top of the plot. UL¬ 
GRB s (black) show pronounced X-ray emission beyond 10 48 erg s -1 for hours 
following the burst, but are still an order of magnitude shorter than the longest 
events (the candidate TDEs with durations of 10 5 - 10 6 s. It is also interesting to 
note that there is one other stellar transient that can achieve similar late time lu¬ 
minosities, the case of the super luminous supernovae SCP 06F6 JLevan et al.| 
|20T3j 


3.4. Multiwavelength afterglow propertjes 
3.4.1. Optical afterglows 

The afterglow properties of the ultra-long bursts form a var¬ 
ied set of properties, but are frequently rather different than 
those observed in the majority of long-GRBs. In the case of 
GRB 101225A, the early optical-UV light curve is exception¬ 
ally blue, and cannot be well fit by fireball models that work 
well in other bursts ( [Thone et aL[|2011[|Campana et al.||2011 ). 
In addition to this it shows marked chromatic evolution from 
blue to red over the course of the next few days, again dif¬ 
ferent from the largely achromatic evolution of afterglow light 
curves. This behaviour is reminiscent of the evolution of the 
low-luminosity, also very long GRB 060218 (Campana et al. 


|2006| ), that is commonly explained as due to shock breakout 
during a supernova. 

The afterglow of GRB 130925A is highly extinguished by 


dust within its host galaxy, with E(B-V) = 2.24+0.26 (Greiner 


et al. 2014[ ). It can be spectrally well fit with an absorbed bro¬ 


ken power-law, without strong evidence of spectral evolution, 


however the results do not appear consistent with other after¬ 
glows, and indicate the the X-rays and optical light may be orig¬ 
inating from different locations (unlike most GRBs where both 
should arise from the forward shock at late times). Perhaps the 
most striking feature of the optical and IR observations of of 
GRB 130925A is that because of its extreme duration, optical 
and IR observations with GROND were possible while the y- 
ray emission was ongoing. These observations appear to show 
the IR light curve with essentially the same structure as seen 


at higher energies, but with a lag of 300-400 s (Greiner et al. 


|2014| ). While rather different from GRB 101225A, this after¬ 
glow also therefore appears highly atypical. 

The afterglows of GRBs 111209A and 121027A on the other 
hand appear to have little extinction, and are well described as 
power-laws with indicies comparable of those found in most 
long-GRBs, although detailed afterglow modelling has yet to 
be attempted, and it is possible that the luminous X-ray emis¬ 
sion will pose problems for combined models in these cases as 
well. 


3.4.2. X-ray afterglows 

All of the ULGRB s drive long-lived and luminous X-ray 
emission (see Figure[4]for a comparison of ULGRB X-ray light 
curves to other events). At X-ray wavelengths the distinction 
between prompt and afterglow emission is not necessarily clear. 
However, throughout the period before the rapid X-ray decay, 
the ongoing flaring, and apparent hardness-intensity correlation 
(e.g. |Butler &~ Kocevski] |2007| ) are consistent with the X-ray 
emission dominated by the prompt component. In this case the 
X-ray afterglows should begin to dominate only after this de¬ 
cline. At this point the luminosity of the afterglows is mixed. In 
the case of GRB 101225A few X-ray’s were observed after the 
rapid decline ( [Campana et al.||20TT] ). For GRBs 111209A and 
121027A, as with the optical a rather typical afterglow was ob¬ 
served ( [Levan et aLj|2014a| , while in the case of GRB 130925A 
the late time X-rays show strong chromatic evolution. This has 
been interpreted by Evans et al.] |2014| ) as an afterglow which 
is almost absent, but with a major contribution from dust scat¬ 
tered light of the earlier X-ray emission. However, in this case 
one would expect very little higher energy emission, in contrast 
to the detection of an underlying hard component by NuSTAR 


( Bellm et al . j2014| that strongly implies a significant afterglow 
component. However, the standard absorbed single power-law 
that provide a good description of most Swift X-ray observa¬ 
tions fail to provide a good description of the data. These mod¬ 


els can work through the NuSTAR regime (e.g. Kouveliotou 


et al.||201~3] ), or can include additional hard power-law compo 


nents as sometimes seen by Fermi-hXI (Zhang et al., 2011). In 
the case of GRB 130925A the combined spectrum actually re¬ 
quires some additional softer component (absorption, thermal 
emission), whose variation compared to the power-law is prob¬ 
ably responsible for the apparent spectral changes seen in the 
Swift-X RT observations. This interpretation is also supported 
by the multi wavelength data obtained by Piro et al. ( 2014| ), who 
also find a good fit with a thermal + power law model while also 
providing a fuller modelling of the afterglow. Although the in¬ 
tensive work on the afterglow of GRB 130925A offers rather 
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disparate interpretations, a common theme of these models is 
that the X-ray afterglow itself is relatively weak, and consis¬ 
tent with being driven into a low density environment. Such 
detailed observations are not available for other events, but the 
absence of the X-ray afterglow in GRB 101225A would also 
support such a model. A consistently low density surround¬ 
ing ULGRBs would provide potentially valuable constraints on 
their progenitors. 


3.5. Host galaxies 

The host galaxies of the ULGRB population are all star 
forming galaxies, with detected emission lines. In the case of 
the first two examples (GRBs 101225A and GRB 111209A) 
the host galaxies were apparently extremely compact and of 
low luminosity (with Mg « -16.2 and -17.5 for 101225A and 
111209A respectively (Levan et al.| 2014a)). These properties 
were consistent with an origin in blue compact dwarfs, with 
likely low metallicity, and such galaxies are rare, even amongst 
the GRB host population. In turn this was taken as some ev¬ 
idence for a distinction between the ULGRBs and the normal 
LGRB population. However, the hosts of GRB 121027A and 
GRB 130925A are rather brighter systems. 

A further important diagnostic is the location of the events 
within the host galaxies. Long-GRBs are highly concentrated 
on their host light (much more so than normal core collapse 
SNe, |Fruchter et al.| ( |2006| ); |Svensson et ak] ( |2010| )), while short- 
GRBs are highly scattered at large offsets ( |Fong et aT. 2010 


2013| ). Again, the first two examples of ULGRBs showed an 


apparent origin consistent with the nuclei of their hosts, per¬ 
haps indicative of an origin with the black hole that may reside 
there. However, these small galaxies also contain a significant 
fraction of their total light within the brightest pixel of their 
hosts (since they are compact), and so the probability of asso¬ 
ciation with this location should they follow a distribution sim¬ 
ilar to GRBs is significant. Indeed, the case of GRB 130925A 
clearly shows that not all events are nuclear, with HST imag¬ 


ing clearly demonstrating an non-nuclear origin (Tanvir et al. 


2013b| ). However, in this case the galaxy also shows signs of 


recent merger activity, and so may harbour more than one mas¬ 
sive black hole. 


deep, late-time multicolour observations were taken in each 
case ( [Thone et al.[ |2011[ |Levan et ah] |2014a| ). In both cases 
there is spectral evolution, with a change from blue —> red. 
This could well be the hallmark of supemovae emission. GRB 
afterglows are typically described as power laws (modified by 
host extinction, or the unusual possibility of a spectral break in 
the range of interest), in contrast SNe exhibit strong metal line 
blanketing short ward of ~ 3000A in the rest-frame. Hence the 
emergence of an associated SNe should provide a measurable 
reddening, and may well reverse the decay at rest-frame optical 
wavelengths where the SNe peaks. 

However, while the first order properties seen in GRB 101225A 
and GRB 111209A appear to support this model, they differ in 
the details. In particular, the counterparts remain too blue to 
be well described by a supernovae similar to the SN 1998bw 
models that provide such a good fit to most LGRBs (see Fig¬ 
ures [5 ]) and spectroscopy taken of GRB 111209A at the time 
of expected SNe peak fails to show any obvious SNe features 
(Figure [6]). 



1Q14 1Q 15 
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3.6. Supernova searches 

Three out of the four well studied ULGRBs lie at redshift 
where supernovae could be visible should they be similar to 
those seen in normal LGRBs (Li et al., 2014| ). Notionally, this 
should be most straightforward in the case of GRB 130925A, 
whose redshift (z = 0.35) is almost identical to the case of GRB 
130427A, whose SNe was readily identified ( |Xu et al.[ |2013| 
|Levan et"aL||2014a| . However, the burst afterglow was appar¬ 
ently highly extinguished ( [Greiner et al.||2014| ), and only seen 
in the near-IR, where detailed follow-up was more difficult. To 
date there is no clear evidence for an SNe in this case. 

This leaves the cases of GRB 101225A and 111209A. Both 
of these events show strong UV emission, and little evidence 
for significant host extinction. Despite their higher redshifts 
they are promising routes to searching for associated SNe, as 


Figure 5: The evolution of the spectral energy distribution of GRB 101225A 
and GRB 111209A (from Levan et al. 2014). Both can be seen to evolve 
from blue to red, with the initial SED of GRB 101225A rising into the blue 
(in contrast to GRB afterglows that scale typical as v _1 ). In broad terms this 
evolution is what is expected from emerging SNe. However, the counterparts 
remain too blue to be well fit with simple SN templates similar to those seen in 
normal GRBs. 


3.7. Progenitor models 

There are various models posited for the origin of ULGRBs, 
each have distinct advantages and constraints when describing 
their observed properties. To date there is no smoking-gun ob¬ 
servation of the progenitor of an ULGRB, and so each of these 
models remain plausible. 

The observed duration of a GRB (T y ) is dependent on the 
time that the central engine is active (T e ), minus the time that 
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GRB 

^90 (s) 

Tburst (s) 

^90, burst (s) 

z 

^iso (erg) 

Host mag 

References 






101225A 

7000 

106659 

21800 

0.85 

5.2 x 10 52 

-16.2 (g) 

Levan et al. 

(2014a 

); 

Thone et al. 

C 

>onl 

111209A 

13000 

63533 

14300 

0.677 

7.0 x 10 52 

-17.5 (r) 

Levan et al. 

(2014a 


Gendre et a] 

.( 

2013) 

121027A 

6000 

35399 

10700 

1.773 

8.8 x 10 52 

N/A 

Levan et al. 

2014a 

); 

Starling ( 

2015 


130925A 

5000 

11614 

2700* 

0.35 

7.0 x 10 51 

-18.4 (g) 

Evans et al. 

2014): 


3 

(201 ; 



Table 1: The Swift sample of ultra-long GRBs with detailed X-ray, optical and IR follow-up. Tgo is based approximately on the duration over which y -ray emission 
was detectable to BAT, but is different from that reported by standard analysis due to the continuing contribution over subsequent orbits. Tburst is taken from |Zhang| 
|et al. 1(2014) , and Tgg^urst is approximately estimated based on the cumulative flux density with Tburst (light cures following |Evans et ak|(2010|f2009) . The redshifts 
given are from the literature, while the isotropic energy releases ( E \ so ) are based on the integrated 10 keV flux density, and do not include any correction to a wider 
band. *In the case of GRB 130925A the brightest emission was missed by BAT, and so the integration yields a smaller Tgo^urst- The host absolute magnitudes are 
given simply by M = m - 5 log(d/10pc)) + 2.5 log(l + z), and are not further k-corrected. The brackets refer to the band of measurement. 



Figure 6: HST grism spectroscopy of the long-GRB 130427A (Levan et al.| 
|2014b) and the ultra-long GRB 111209A (Levan et ak]|2014a) . The flux levels 
have been scaled to match for comparison of the overall spectral shape, and 
have not been corrected for either afterglow or host galaxy contribution. The 
case of GRB 130427A clearly shows features consistent with a SN 1998bw SNe 
(in particular a broad bump at around 5000A). In contrast, no such features 
are seen in the lower signal to noise spectrum of GRB 111209A. However, 
despite the signal to noise it is apparent that the overall spectral shape in GRB 
111209A is much flatter (in Ff) than expected for an SNe. This could be due to 
a larger contribution from the host galaxy, because the SNe type is different (for 
example the flatness in this spectral regime may be due to enhanced emission 
around 6500Afrom a broad Her in a SN II), or because there is no SNe in this 
ULGRB. 


the jet takes to penetrate whatever dense medium surrounds the 
central engine ( T b ), i.e. T y -T e -T b (|Bromberg et al.||2012). 
Ultra-long bursts then require pro-longed activity of the cen¬ 
tral engine. For most GRBs it seems likely that the progenitors 
are very compact, leading to small breakout times. Fine tun¬ 
ing issues would be problematic in creating short periods of 
y-ray emission if the time to breakout was much longer (i.e. the 
scenario where T e « T h » T y would require fine tuning). For 
ULGRBs such constraints are not present, and it is possible, but 
by no means required, that much longer breakout times (corre¬ 


sponding to much larger radii progenitors) can create ULGRB s. 


3.7.1. Collapsars 

Collapsars are appealing as the progenitors of ULGRB s since 
they are known to produce long duration GRBs and so can 
clearly create events with the necessary energy, spectra etc. It is 
less clear how they may create such long duration events, and 
what may distinguish the longest duration events from “nor¬ 
mal” LGRBs. In most GRBs the long-lived engine is thought 
to be powered by fail-back from material expelled in the burst, 
but with insufficient velocity to escape the newly formed black 
hole (e.g. |Zhang et al.||2006l|Fryer et al.[|2Q07| ), or alternatively 
from material that fragments during the collapse process ( |King| 
et al. 2005| . In the former case one would expect the infall rate 
to reduce with tim e (naively the fallback accre tion rate scales 
asMoc r 5/3 (e.g. Fryer 2006 Fryer et al. 2007 1). In the latter 
case one might expect flares, rather than continuous emission. 
This doesn’t easily produce the observed light curve properties 
of ULGRBs. 

However, the observed y-ray s are not just the product of the 
central engine, but (at least in the case of LGRBs) the inter¬ 
action of the shocks that are produced in the outflow from this 
engine. These shocks decelerate in the medium surrounding the 
burst, and then shock with it to create the afterglow. It has been 
suggested in the case of GRB 130925A (and possibly applica¬ 
ble to other ULGRBs) that a low density medium surrounding 
the GRB would result in a significantly longer deceleration time 
(Evans et al. - ;2014[ ), with internal shocks between ejected shells 
continuing to produce y-rays over this period. The challenge for 
this scenario is that the low density medium would be expected 
to drive only a weak afterglow. In the case of GRB 130925A 
this is side-stepped since it appears to be in a dusty environ¬ 
ment, and the unusual afterglow that softens with time may be 
explained by a dust echo of the prompt emission (Evans et al 


2014). It may also not represent a problem for GRB 101225A, 
where there is essentially no X-ray emission beyond the early 
rapid decay (and so potentially no afterglow). However, in the 
case of GRB 111209A and GRB 121027A, these is little evi¬ 
dence for significant dust extinction, and the afterglows appear 
much more typical (Le van et al., 2014a ). 

Finally, the absence of obvious SNe signatures in the af¬ 
terglow of GRB 111209A would appear to conflict with this 
model. With two exceptions, that may belong to a short-population 
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of bursts (e.g. |Fynbo et al.[|2006[ ), LGRBs are known to be as¬ 
sociated with high velocity type Ic supernovae. Indeed, recent 
work shows these to be remarkably uniform, and may even en¬ 
able their use as standard candles e et~aE| 2014 ). However, 
in the case of GRB 111209A there is no such signature, de¬ 
spite deep HST observations that should have readily detected 
it (|Levan et al.’ |2014a| ) (see Figure [5]). This implies that this 
burst was not a typical LGRB, and hence its progenitor may 
well have not been a typical collapsar. 


3.7.2. Collapse of giants or supergiants 

The collapse of giant stars offers an alternative scenario to 
explain ULGRBs ( [Gendre et al.| |2013[ |Levan et aLj |2014a ). 
These systems are known to create type IIP SNe based on the 


direct detection of their progenitors (e.g. Smartt et al. 2009), 


but their ability to create high-energy transients remains un¬ 
clear. While collapsar progenitors typically have radii of only 
a few Rq , supergiants have radii of an AU or greater. Their 
collapse times are therefore orders of magnitude larger than in 
collapsars, and given their lower densities their M rates are also 
smaller. The fuelling of the nascent black hole would then nat¬ 
urally create a burst of extremely long duration. However, it 
is less clear if the collapse would necessarily create the neces¬ 
sary accretion disc geometry to form a GRB-engine, or if the 
Eddington limit should start to apply during the prolonged col¬ 
lapse, essentially halting the infall and stalling the GRB. How¬ 
ever, numerical simulations (Nakauc hi et al. 2013| ) do appear 
able to drive transients with the requisite durations from giant 
star collapses. 

GRB 130925A offered an ideal opportunity to test these 
models, which were largely developed as a result of observa¬ 
tions of the earlier bursts. Unfortunately optical extinction from 
the host makes a direct search for SNe II impossible, however, 
harder X-ray observations do not suffer from the same problem 
and have been exploited to this end. While Evans et al. |2014| 
interpret the evolution of the Swift-XRT data as due to an evolv¬ 
ing dust echo, the combined Swift , XMM-Newton and NuSTAR 
observations are well fit with a power-law plus a thermal com¬ 
ponent. |Bellm et ah] p014| assume this thermal component to 
originate in a hypercritically accreting disc, whose luminosity 
is consistent with being fed by fallback from a supergiant star 
(although the range of allow fallback rates is large enough to en¬ 
compass many models ( Wo ng et ah} 2014| )) . Piro et al. ( |2014 ) 
advocate an alternative interpretation in which the thermal com¬ 
ponent arises from cocoon emission surrounding the relativistic 
jet, and importantly note the requirement of a low ambient den¬ 
sity, consistent with one sculpted by a blue supergiant star at 
low metallicity. 

While it remains unclear if this model will survive further 
scrutiny (for example the the host region appears to be at high 
metallicity ( [Schady et al.|[2015l ), and the robustness and unique¬ 
ness of the thermal solution is unclear (Ev ans et ak]|2QI4| )), fu¬ 
ture observations both in the X-ray and at optical wavelengths 
should provide direct tests of giant and supergiant models. 


3.7.3. Tidal disruption flares 

An alternative model for the origin of ULGRBs is that they 
arise from tidal disruption events. The appeal of this model is 
two-fold. Firstly, a tidal model has been invoked with some suc¬ 
cess for the very longest transients (see section]?]), and secondly, 
the first two events with detailed follow-up (GRB 101225A and 
GRB 111209A) appeared to lie in locations consistent with the 
nuclei of their hosts in deep imaging ( Levan et al] 12014a[ ). 

However, classical versions of the TDF models in which 
main sequence stars are disrupted by supermassive black holes 
rapidly run into problems. Importantly, the timescale of UL¬ 
GRBs of ~ 10 4 is significantly shorter than the orbital period 
at the tidal radius of a 10 6 M© black hole and main sequence 
star. Secondly, the locations of other events, in particular GRB 
130925A are clearly non-nuclear, and so would require unusual 
merging systems to create (although it should be noted that in 
the case of GRB 130925A this is a distinct possibility). 

A more appealing prospect is that these events do not repre¬ 
sent normal TDFs, but rather variants in which lower mass, in¬ 


termediate mass black holes tidally shred white dwarfs (Krolik 


|& Piran[|2011[[Levan et al.[|2014a[ [MacLeod et al.[|2014| ). At 

these low masses the tidal radius for a white-dwarf can lie out¬ 
side the Schwarschild radius, something impossible for more 
massive black holes. The compact orbit naturally matches the 
timescale observed in ULGRBs, and, unlike SMBH systems, 
it might be expected that IMBH could arise in young stellar 
clusters, where they may also power Ultra- or Hyper-luminous 
X-ray sources, or in the cores of globular clusters, all of which 
would place them well away from the nuclear regions of more 
massive galaxies, but perhaps in the nuclei of low mass systems, 
such as the hosts of GRBs 101225A and 111209A. 

3.8. Future prospects 

The combination of extremely long durations, with frequently 
unusual and unexpected multi wavelength properties offers a 
strong guide that the ULGRB population does represent a gen¬ 
uinely new class of burst. The crucial question underlying stud¬ 
ies of the ULGRB population is then the unveiling of the pro¬ 
genitors. This diagnostic would both directly answer the ques¬ 
tion of whether ULGRBs are genuinely a new population, or 
the long tail of the ULGRB distribution, and identify how they 
might be used as cosmological probes. If they were somehow 
related to SN Ic as normal LGRBs are then it would clearly be 
important to understand what features in the core collapse give 
rise to such long-lived and luminous central engines. Models 
that rely to date on fallback to create plateaus or on fragmenta¬ 
tion of the supernovae ( [King et aL]|2005| ), or in the disc ( |Pema| 


et al., 2006) would seem to fail to provide the necessary late 


time accretion rates, and so our understanding of core-collapse 
in GRB environments would seem in need of revision. 

Alternatively, should they arise from the collapse of giant 
stars then they provide a demonstration that engines can be ac¬ 
tive in a far greater range of supernovae than currently appre¬ 
ciated. The ability of these engines to power GRB-like events 
in massive stars that retain their hydrogen envelopes may make 
ULGRBs relatively common through the Universe, while their 
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observed rarity is largely due to selection effects (Levan et al., 
2014a). Indeed, the ability for engines to be active in relatively 
extended stars may bode well for their use as probes of popula¬ 
tion III stars ( [Piro et al^|2014| ). 

Finally, should they in fact arise from tidal disruptions in¬ 
volving intermediate mass black holes then they would pro¬ 
vide one of the few robust routes to the identification of such 
systems, and in turn enable uses as probes of the locations of 
IMBHs in dwarf galaxies and globular clusters, as a route to 
hone the M - cr relation, and as an input to models of galaxy 
formation and evolution. 

It is interesting to note that all of the above systems seem 
plausible for the origin of ULGRBs, both in the sense that the 
rates of the relevant progenitor systems mean they should oc¬ 
cur, and because more detailed simulations manage to create 
transients with the relevant durations in each case ( |Nakauchi| 
et al.||2013}|MacLeod et al.||2014| ). It may then be the case that 
more than one of the above scenarios provide an explanation for 
ULGRBs (or at least long duration high energy transients more 
generally), or that there are alternative routes of identifying the 
other channels. 

The actual task of making the distinction between these 
models should be relatively straightforward with current tech¬ 
nology. SNe in GRBs can be observed from the ground to z ~ 1, 
and potentially beyond with HST. This hypothesis is already 
stretched by the existing observations, but can be rigorously 
tested with a modest set of new observations. The SNe asso¬ 
ciated with the collapse of giants should be hydrogen rich, but 
could be an order of magnitude fainter than those seen in GRBs. 
This makes their detection more difficult, but the presence of 
strong, variable H a emission should make them traceable out 
to similar redshifts with current IR spectrographs, or again with 
HST observations. Finally, in white dwarf disruptions there may 
be no SNe, or an unusual type I event (different from the type 
Ic associated with LGRBs ( [MacLeod et al.[|2014| )). It is inter¬ 
esting to note that the origins of both LGRBs, and recently of 
SGRBs have been pinned down by studies of late, red bumps 
in their optical and IR light curves (e.g. Hjorth & Bloom, 2012; 


Cano| |2013| |Tanvir et ak| |2013a[ |Berger et al.| 2013| ). These 

studies are likely to be equally diagnostic for ULGRBs, and 
should be possible in the next few years. 


4. The longest events, beyond 10 4 5 s 

As well as the population of ULGRBs, that have durations 

of T 90 ~ 10 4 * * * s, there is an additional population of newly 
recognised y-ray transients that are detectable for days after 

the initial trigger (see Figures [T] and [7]). These have been de¬ 

tected initially as GRBs, but also via the BAT transient moni¬ 
tor programme ( [Krimm et al.[ |2013| that reconstructed the im¬ 
age plane for BAT and searches for transients on even longer 
timescales. While some very long y-ray outbursts arise from 

either Galactic binaries, or well known AGN, this new popula¬ 
tion do not. Instead they occur at cosmological distances, reach 
peak luminosities in excess of 10 48 ergs s -1 (the Eddington limit 

of a 10 10 M© black hole), but are found in galaxies showing lit¬ 
tle or no evidence for AGN activity. 



Figure 7: Cartoon phase space for various classes of high energy transient of 
differing durations, modified from |Levan et al. The axes shown an 

X-ray and optical luminosity (expressed as an absolute magnitude), while the 
background points arise from a cross correlation of the 2-XMM catalog jWat-| 
|son et al.] 2009 > with SDSS ( jAbazajian et al.| |2009f and a sample of local 
galaxies observed with XMM-Newton lHo] 20094. Each population (LGRB), 
(ULGRB), (rTDF) is shown at a typical peak time (100s for GRBs, 10 4 s for 
ULGRBs and 10 5 s for rTDFs). The arrows indicate the evolution of the events 
to 10 6 s. There is significant diversity in the LGRB population, and this may 
well be represented as well in the overall populations of ULGRBs and rTDFs. 
Therefore both the locations of the populations, and their evolution should be 
viewed as approximate. Nonetheless, it is clear the events do occupy quite dif¬ 
ferent regions of phase space, and diagnostics such as this may be powerful 
discriminants between differing origins for transients, in particular the TDFs. 


These events are generally ascribed as relativistic variants 
on tidal disruption flares, in which some small fraction of the 
ejecta is emitted at relativistic speed. This model is compelling, 
but alternatives have been suggested. Perhaps unsurprisingly 
these suggestions are very similar to those discussed for the 
origin of ULGRBs. In particular this includes both WD-IMBH 
disruptions, and long lived accretion in very massive, extended 
stars. 


4.1. The sample: 

Swift J1644+5734, Swift J2058+0516, Swift J1112-8238 


The best studied candidates TDF is Swift J1644+5734 (Leva i 


et al.l|2Ql7||Bloom et al.||20lTl|Burrows et al.l|2011[|Zauderer 


et al. , |2011| ), which was detected via a GRB trigger and had sub 
stantial additional follow-up. It was initially thought to be likely 
due to a soft, fast X-ray transient with the Galaxy, and indeed 
similarities were noted with the Galactic system IGR J16479- 
4514 ( [Kennea et al.[ |20TT| ). However, redshift measurements 
from Gemini and the GTC ( [Levan et al.[ |2011| ) demonstrated 
an extragalactic origin at z = 0.35, and highlighted the extreme 
luminosity of the event. In total Swift J1644+5734 triggered the 
Swift-BAI on four occasions over 48 hours, and went on to be 
detectable by BAT for over a week (Levan et al. [201 If Burrows 
|etaLl|2QTTT ). 


The second example located is that of Swift J2058+0516 
( Cenko et al.[|2012) , which was detected by the Swift transient 
monitor only 2 months after the detection of J1644+5734. 
Follow-up observations again revealed a cosmological source, 


9 





































































this time at z = 1.18, based purely on absorption lines in Keck 
spectroscopy (no emission lines were seen in this, or in X- 


shooter observations covering the common host lines. Pasham 
letH] ( [20l5] )). 

The final example to date, Swift J1112-8238 was not located 
in real time, but in an archival search within the BAT transient 
monitor archive ( |Krimm et al.[|2013l|Brown et al.[|2015[ ). This 
means that real-time follow-up was limited, but later observa¬ 
tions reveal an X-ray and optical transient coincident with the 
nucleus of a faint (and hence low signal to noise) galaxy, likely 
at z = 0.89 base o n a s ingle emission line identified as the 0[n] 
3727A doublet (Brown et al. 2015). 


4.2. X-ray counterparts 

The X-ray counterparts of these transients are all extremely 
luminous and long-lived. They have peak luminosities around 
10 48 ergs s -1 , and remain brighter than 10 46 ergs s -1 for 50- 
100 days. Indeed, all three have remarkably similar luminosi¬ 
ties at around the 50 day mark. It is difficult to define the 
power-law slope precisely in these cases, since the zero time 
of the burst is unclear. For example, in the best studied case 
of Swift J1644+5734, there were multiple triggers of BAT over 
several days, but detailed searches revealed the source was ac¬ 
tually detected at least 4-days before the first of these ( Burrows 
|etaO 20 Tl). Additionally, the light curve shows pronounced 
dips, that have been claimed so show some evidence for peri¬ 
odicity (Saxton et al., 2012| ). This makes fitting a non-trivial 
task, although a simple fit to the late time light curve of Swift 
J1644+5734^from 10 6 to 4x 10 7 syields an index of t~ 165 . The 
case of Swift J2058+0516 appears to be rather steeper than this, 
with a decay of around t~ 2,2 (Cenko et al., 2012). Finally, for 
Swift J1112-8238, the sparse data provid es a limited lever arm 
but is best fit with a decay around t~ 1A (Brown et al. 2015). 


These decay rates can in principle be compared to the expecta¬ 
tions of tidal disruption models, although recent detailed calcu¬ 
lations suggest that the simple approaches previously taken may 
be very different from what will actually be observed, although 
the t 5/ 3 and t~ 22 of Swift J 1644+5734 and Swift J2058+0516 are 
remarkably close to certain scenarios (Guillochon & Ramirez- 
|Rmzl|20l3) . 

In addition to this broad decay there are two additional im¬ 
portant features in the X-ray light curves. The first is rapid 
variability occurring at early and late times, often showing a 
doubling of the flux on timescales of only a few hundred sec¬ 
onds, even hundreds of days after the burst (L evan et al., |2Q11 1 
Pasha nTet al.[ |2015| ). This behaviour is similar to that seen in 
blazers, and so is likely a hallmark of a relativistic jet. The 
second striking feature observed in Swift J1644+5734, Swift 
J2058+0516 and possibly Swift J1112-8238 is a rapid shut off 
of the X-rays, a year or more after the outburst, in the case 
of Swift J1644+5734, the counterpart appeared to switch off 
within a few days, interpreted _as a power-law slope it would 
have been steeper than r 70 (Levan et al., 2015), see Figure [8] 


This combination of X-ray properties (luminosity, variabil¬ 
ity, shut-off) are remarkable, and are not mirrored in any other 
class of source. 

4.3. Optical and IR observations 

Optical observations are rather more limited than those in 
the X-ray for these sources, especially in the case of Swift J1644+5734, 
which had nearly daily monitoring at X-ray wavelengths in the 
months (and even years) since its discovery. However, optical 
counterparts have been discovered to all three of the candidate 
events. These vary in their apparent properties. In the case of 
Swift J1644+5734 the source was highly reddened with Ay ~ 6, 
not surprising if it lies in, or close to the host galaxy nucleus 
( [Levan et al.[|201l||Bloom et al.[|20lH|Burrows et al.||2011| ). 

In the case of Swift J1112-8238 only a handful of detections 
were obtained at the time, since the source was not identified as 
of interest in real time. Only for Swift J2058+0516 do we have 
a good, and well sampled optical/IR light curve from early to 
late times after the outburst. This shows clear evidence for a 
hot, thermal component that cools at near constant radius in the 
months and years since the outburst. This would appear to be 
distinct from the non-thermal X-ray emission, perhaps indica¬ 
tive of a different location for its production. However, it is 
interesting to note that it does apparently switch-off at the same 
time as the X-ray emission, perhaps indicating a similar origin, 
despite the rather disparate spectra ( [Pasham et al.||2015| . 

Optical observations provide the strongest possible constraints 
on the positions of the sources within their host galaxies. In the 
cases of Swift J1644+5734 and Swift J2058+0516 this is possi¬ 
ble thanks to early and late time observations with HST ( |Levan| 

|et al.| |2011| |Pasham et al. 2015| ), and is shown graphically in 
Figure [9] In these cases the sources are shown to lie within 150 
and 400 pc of the nuclei of their host galaxies in projection, and 
in each case are consistent with a nuclear origin. In the case of 
Swift J1112-8238 only ground based images are available at the 
time of writing. These show the source to be associated with the 
centroid of its host galaxy, but the low signal to noise precludes 
strong conclusions from this. 

4.4. Radio counterparts 

One of the most striking features of Swift J1644+5734 is 
the presence of a bright and extremely long-lived radio counter¬ 


part (ZaudereretaL 2011[|Berger et al.||2012||Zauderer et al. 


2 taken from, http://www.swift.ac.uk/xrt_curves/450158 


2013). This rises on a timescale of 200 days after the burst, 
continuing to brighten well beyond the X-ray peak. Indeed, 
it is still visible beyond the point at which the X-rays appear 
to shut-off (Zauder er et al. 2013| ). The radio counterpart pro¬ 
vides direct evidence for relativistic motion, although with only 
a modest Lorentz factor T ~ 2 Zaude rer"et al.| ( [20TT] ). Should 
the jet associated with these outbursts be similar in behaviour to 
GRB jets this would imply that it is relatively wide. Limited ra¬ 
dio follow-up for Swift J2058+0516 would seem to match these 
expectations, although much more sparsely sampled. 

A striking feature of these radio luminosities is that they 
should be visible to next generation radio surveys to z > 6. 
Combined with the longevity of the luminous radio emission 
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(several years at z ~ 6) it may well be that these sources will 
ultimately be best located in the radio regime. 



Figure 8: The X-ray light curve of Swift J1644+57 on a linear scale, highlight¬ 
ing the rapidity of the switch off. The red line prior to the switch off is a t~ 5 ^ 
decay, while the extension of this (two orders of magnitude fainter) is shown 
as the dashed line beneath. At present it is not clear if the level seen after the 
drop-off is some continuation of a decay, or a quiescent level in the host galaxy 
(e.g. from an underlying AGN). 


4.5. Relativistic tidal disruption flares 

The concept of a tidal disruption flare has long been dis¬ 
cussed, and is an inevitable consequence of the presence of su- 
permassive black holes in galactic nuclei. The disruption occurs 
when a star approaches a black hole sufficiently closely that the 
gravitational pull on its outer layers from the black hole exceeds 
that from the star itself. At this point the star can disrupt either 
completely, or partly (e.g. by stripping of the outer layers of 
the star). In a complete disruption half of the star is placed onto 
eccentric bound orbits and will return to form a disc and ulti¬ 
mately accrete onto the black hole, the other half is unbound. 
The accreting material is shock heated to high temperatures 


peaking in the extreme-UV or soft X-ray regimes (e.g. Rees 


|1988| ). Numerous candidates for such events have been uncov¬ 
ered (e.g. Gezari et al. 2003 [Komossa et ah) 2004[ |Halpern| 
et al.| 2004 |Esquej et al. 2008| " Saxton et al. " |2012 ), although 


it has been difficult to rule out alternative explanations, such as 
unusual nuclear supernovae in many cases. 

However, more recently it was suggested that in addition to 
this thermal emission a small fraction of the material may be 
expelled at relativistic velocities ( Giannios & Metzger) 2011} 
| van Velzen et al.| |2011| ). Such material would create a new 
route to identifying TDFs, and was originally considered in the 
context of next generation radio surveys, which could uncover 
the events at late times when the blast wave was approximately 
spherical. 

These extremely long transients provide a natural match to 
the expectations of this model. There is direct evidence from 
the radio observations for relativistic motion, while the super- 
Eddington accretion rates observed at early times might suggest 


an even narrower beam (the early luminosity exceeds the Ed¬ 
dington limit for the black holes within these host galaxies by 
a factor of > 10 4 ). All of the events are consistent with the nu¬ 
clei of relatively low mass galaxies, and with black holes with 
M B h < 10 8 M© (see below), consistent with the disruption of 
main sequence stars occurring outside the event horizon. Given 
the theoretical prediction of these events prior to their detection, 
their long-lived, high luminous emission, and their locations 
within host galaxies this model provides an extremely good fit 
to the observations, and is adopted as the consensus view for 
the origin of these outbursts. However, it is not a unique model, 
and other alternative has been put forward as well. 


4.5.1. Black Hole Masses 

If these events are interpreted as TDFs then various routes 
have been proposed to measure the black hole masses within 
them. The most straightforward route (though potentially with 
the largest scatter) is to utilise the well known relationship be¬ 
tween bulge-mass and black hole-mass in galaxies (e.g. |Haring] 
& Rix, 2004| ). The total stellar mass of the galaxies can be es¬ 


timated from rest-frame IR observations, and then used to infer 
a black hole mass. In each case, the masses of the black holes 
obtained via this route are M BH < 2 x 10 7 M 0 . 

Alternatively, one can use timing arguments to infer the 
mass of the black hole. The rapid variability of Swift J1644+5734 
implies that any temporally connected source by at must be 
smaller than the Schwarschild radius of a 8xl0 6 M© BH (Bloom 


et al. 201 1| ), providing a relatively robust upper limit to the 


black hole mass. Alternatively, a rather more complex approach 
can utilise known relations between radio and X-ray luminos¬ 
ity, that apparently links stellar and supermassive black holes. 
This route infers are rather small intermediate mass black hole 
lo gM BH = 5.5 ± 1.1 M© |Miller & Gultekin| ( [2QTTT ). 

Finally, it has been suggested that the rapid switch off seen 
in both Swift J1644+5734 ( [Levan & Tanvir 2012 |Zauderer| 
et al.||2013] ) and Swift J2058+0516 ( Pasham et al. 2015 ) might 


be due to the source reaching its Eddington limit. This has 
some appeal, since at this point it would no longer be neces¬ 
sary for the source to be beamed in order to drive its luminosity 
against inflating material. In this case, the luminosity at the 
break would be indicative of the the Eddington limit for that 
mass of black hole, providing a route to measuring the mass. 
A conservative approach would be to assume that the Edding¬ 
ton limit lies somewhere between the pre- and post-drop lumi¬ 
nosity, although the drops are so steep that this only provides 
a black hole mass limit in the range 10 4 < M bh < 10 6 M© 
( [Pasham et al.||2015| . 

All of these mass estimates contain considerably uncertainty, 
and rely either on the extrapolation of known relations signifi¬ 
cantly beyond the regions in which they have been tested (e.g. 
the M\> u i ge -M BH relation, or the variability fundamental plane), 
or on hitherto untested relations (i.e. that the jet switch-off oc¬ 
curs at the Eddington limit). However, it is encouraging that 
all of these routes yield masses of M BH < 10 8 M© for all three 
events, meaning that they are consistent with the expectations of 
main sequence stars being shredded by black holes, comparable 
in mass to those that we know exist within the local Universe. 
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GRB 

Trigger 

z 

peak L x (erg) 

Mrh (M 0 ) 

References 




Swift J1644+5734 

Burst 

0.35 

3 x 10 4S 

10 5 < Mpu < 10 7 

Levan et al. ( 

2011); 

Burrows et al. 

(2011) 

Swift J2058+0516 

Monitor 

1.19 

7 x 10 47 

10 4 < M bh < 3 x 10 7 

Cenko et al. 

2012); 

Pasham et al. 

(2015) 

Swift Jl 112-8238 

Monitor 

0.89 

7 x 10 46 * 

~ 2 x 10 7 

Brown et al. 

(2015) 


Table 2: The Swift sample of transients with duration > 10 5 s, the candidate relativistic tidal disruption flares. * indicates that X-ray observations started very late. 
The extrapolation back to early times using the BAT transient monitor would have resulting in a peak luminosity approximately a factor of 10 brighter. The black 
hole mass estimates are based on a bulge luminosity to black hole mass relation a at the upper end, and on scaling arguments relating to the variability or rapid 
shut-off at the low mass end. 


4.6. White dwarf disruptions 

The early light curve of Swift J1644+5734 is punctuated 
by a series of extremely rapid flares, during which the lumi¬ 
nosity jumps by orders of magnitude on a timescale of only 
a few hundred seconds. These flares, repeat on timescales of 
around 30,000 s (although it should be stressed that they are 
not strictly periodic) and have a marked similarity flare to flare. 
These timescales are perhaps surprising for a main sequence 
disruption, since they are much shorter than the orbital period 
of the disrupting star, and therefore consecutive passages of the 
star (with further disruption) cannot be explanation for the ori¬ 
gin of the flares. As an alternative, it has been proposed that 
Swift J1644+5734 is due to a white-dwarf disrupted by a IMBH 
( [Krolik & Piran[|2011 ), an interestingly similar model to those 
subsequently proposed by some authors to explain the ULGRBs 
( [MacLeod et al.|[2Ql4| ). The appeal of this model lies in the nat¬ 
ural interpretation of the flare structure. However, interpreting 
this structure is challenging as it likely based on the beamed 
X-ray emission, and so its direct correlation to the M onto the 
black hole is less clear. Furthermore, while most models would 
predict that there is relatively smooth return of material from 
the disrupted star back to the black hole, this may well not be 
the case at the point of return of the most bound debris (i.e. the 
material that first accretes), which could fragment or accrete in 
sporadic bursts, even from a main sequence star. 

4. 7. Massive star collapse 

Swift J1644+5734 was first discovered as a GRB, and its 
long-lived emission was subsequently interpreted as being re¬ 
lated to a massive star collapse (Quataert & Kasen, 2012). More 
detailed models have subsequently been derived by Woosley 
& Heger ( Woosley & Heger, |20l2] ), who demonstrate various 
means by which long-lived transients can arise from collapsar- 
like events, including single massive stars, tidally locked bina¬ 
ries, and pair instability collapse. Interestingly, these models 
predict a rapid switch-off in the in accretion rate, essentially 
when all of the star has been accreted. This may mirror the 
switch-off’s seen in the events to date, and it is striking that 
this prediction was made prior to their discovery (see Figure 1 
of IQuataert & K asen ( 2012)). Indeed, further support that core 
collapse events can reach such high luminosities comes from X- 
ray observations of the super luminous supernovae SCP 06F6, 
that reached L x ~ 10 45 erg s -1 , 150 days after its initial outburst 
( [Levan et al.[|2013| ) 

However, the locations of the bursts are problematic for 
massive star collapse models. Even GRBs, that are highly con¬ 


centrated on their host like only lie coincident with the nucleus 
of their hosts 10-15% of the time ( [Fruchter et al.|12006}|Svens-| 
2010| ). Given this the probability of 3 events having 


son et al. 


such a location is < 0.4%, apparently ruling out the possibility 
of a GRB-like origin at > 3<x. However, it should be noted that 
this conclusion is currently only robust in two cases, and deep 
HST observations of Swzyi Jill 2-8238 are hence crucial to this 
conclusion. 
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Figure 9: Hubble Space Telescope imaging of the host galaxies of Swift 
J1644+57 (Levan et al.||20TT) and Swift J2058+0516 (Pasha m et al.| [20l5> 
in the F160W IR band. The top panel is an early image containing the after¬ 
glow and host galaxy light, the second an image showing just the host galaxy, 
while the third image shows the subtraction of the two epochs. In each case it 
is clear that the sources lie at locations coincident with the nuclei of their hosts. 


4.8. Role as probes 

While the origin of the ULGRBs remains an open question, 
it seems likely (if not yet conclusive) that the longest events 
are due to the tidal disruption of stars by massive black holes. 
These systems are apparently much rarer (in terms of volumet¬ 
ric rate) than “classical”, non-relativistic tidal flares^which are 
expected at rates of ~ 10 -5 yr -1 galaxy -1 (e.g. Rees 1988). 
This is likely because of either the narrowness of the beamed 
emission (so that only a small fraction are seen) or because of 
the rarity of jets (because few systems launch relativistic jets 
in the first instance). Crucial questions therefore arise in the 
role of these events as probes of jet formation, launching an 
ubiquity around massive black holes. The tidal flares are un¬ 
usual because they provide the opportunity to witness accretion 
around a supermassive black hole from its start to finish over 
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the timescale of only a few years, therefore providing new con¬ 
straints about how accretion behaves in a range of regimes. 

Furthermore, since these events may be visible either to X- 
ray or radio searches to much larger distances than non-relativistic 
TDFs, it is possible that they may come to dominate the number 
of systems observed. If a route can be found to translate their 
observations into black holes masses this provides the ability 
to track changes in the black hole mass function through on¬ 
going accretion and heirarchical mergers over most of cosmic 
history. Indeed, it has been suggested that merging black holes 
may raise the tidal disruption rate by 4 orders of magnitude (to 
one per 10 years or so), in which case their discovery might 
provide a route of identifying close merger black hole pairs 
of interest to low frequency gravitational wave detectors (e.g. 
eLISA) (e.g lStone & Loe5H2011 ). 

Finally, the strong particle acceleration in the jets of rela¬ 
tivistic tidal disruption flares are promising sites for the accel¬ 
eration of ultra-high energy cosmic rays (Bloom et al.. 2011 1 
Cenko et al., 2012). Since their are insufficient bright AGN 


within the local horizon through which the rays can propagate, 
and GRBs appear to be insufficient to provide this acceleration 
( |Abbasi et al. 2012| ), it is possible that the tidal flares could 
provide the solution, although further studies of the details of 
shock acceleration, and tidal flare rate are clearly necessary to 
firm this conclusion up (Farrar & Piran , |20T4] ). 


5. Open questions and future prospects 

It is remarkable that highly energetic transient events (E iso > 
10 53 erg) with durations from hours to days apparently exist, but 
have been undetected by a fleet of both pointed and survey y- 
and X-ray telescopes over the past 5 decades. They have evaded 
detection by a combination of rarity (hence they have not been 
found serendipitously in the narrow field of view of X-ray tele¬ 
scopes) and low peak-flux (so they haven’t, in general, been 
identified by rate based y-ray detectors). These events are only 
now uncovered in sufficient numbers to identify new popula¬ 
tions of outburst, thanks to the ability of Swift to both locate 
them, but also rapidly highlight their interest to the community 
in a timely way. However, these new populations are a recent 
development, and our knowledge and understanding is much 
less mature than for the LGRBs and SGRBs, whose origin and 
physics have been hotly discussed for several decades. Future 
observations will continue to build on our knowledge of these 
populations, and address several central questions, including. 

• Are ULGRBs are separate class, spawned from giant stars, 
or WD-IMBH disruptions, or do they form the tail end of 
the LGRB distribution? 

• Can a larger sample of candidate relativistic TDFs con¬ 
firm that they are all related to tidal flares, and not other 
possible progenitors? 

• How large are the true volumetric rates of the systems? 

• What are their roles as probes of both cosmology (pop 
Ill-like ULGRBs, ubiquity of black holes in galaxies) and 


of extreme physics (can their engines explain the origin 
of ultra-high energy cosmic rays? What do their engines 
tell us about jet-formation in AGN?). 

The study of these extreme, but still mysterious transients 
remains an open question in astrophysics, it will be at the fore¬ 
front of work as Swift enters its second decade of operations. 
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